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Abstract 



It is shown that the new data of the DAMA/Nal detector, combined with 
recent calculations of nuclear structure functions, provide the most stringent 
upper bounds on the WIMP-nucleon cross-section both for spin-dependent 
and for scalar interactions. It is proved that, for the first time, data from an 
experiment of direct search for WIMPs allow the exploration of a region of the 
supersymmetric parameter space not yet fully investigated at accelerators. 
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I. INTRODUCTION 



An experiment for direct detection of Weakly Interacting Massive Particle (WIMP), 
which employs a large-mass, low-background Nal(Tl) detector at the Gran Sasso Laboratory 

has recently allowed a substantial improvement of the upper limit on the WIMP-nucleon 
elastic scattering cross-section. A high exposure of 4123 Kg day, a strong reduction of the 
residual internal contaminants and a statistical pulse shape discrimination have made this 
result possible. 

In this letter we examine the impact that this new measurement has in terms of relic 
neutralinos, and more generally in terms of an exploration of the supersymmetric param- 
eter space. First, in Sect. 2 we discuss how upper bounds on (scalar and spin-dependent) 
WIMP-nucleon cross-sections are extracted from the experimental counting rates. Whereas 
the procedure is straightforward for the scalar case, for the spin-dependent case appropri- 
ate specifications on the interaction mechanism and realistic nuclear-physics calculations 
are needed. We show how the new data from the DAMA/Nal-detector provide the most 
stringent upper bounds both to scalar and to spin-dependent WIMP-nucleon cross- sect ions 
(except for a very narrow range of the WIMP mass around 15 GeV). This remarkable result 
derives from the significant improvement in sentivity in the experiment of Ref. [|l| and, for 
the spin-dependent part, also from new theoretical calculations of nuclear effects 0,^. Then, 
in Sect. 3 we proceed to a theoretical analysis of the implications that the new experimental 
data have for relic neutralinos, and we show that, for the first time, data from an experi- 
ment of direct search for WIMPs allow the exploration of some regions of the supersymmetric 
parameter space, whose investigation at accelerators is still under way. 



II. RATES FOR DIRECT DETECTION 



The experimental activity in the direct search for relic particles has been very intense 
in recent years f^^lTl, with quite significant improvements in various technical aspects that 
have drawn the experimental sensitivity nearer and nearer to the one required to match 
the predictions of realistic theoretical models. It is well known that the only safe way of 
discriminating a real signal for WIMPs from the background is to measure a modulation 
effect . Some preliminary measurements |ll2| , p!6| , p!7| have already been performed, but the 
search for modulation will have to be pursued with increased sensitivity in order to reach 
the level required by a number of specific dark matter candidates (such as the neutralino). 
Since the current counting rates are still expected to be dominated by background effects, 
only upper bounds on WIMP-nucleus cross-sections can be derived at present. Some of the 
considerations presented in this section apply to a generic WIMP, however, since we will 
concentrate on the neutralino (x) later on, when we compare experimental bounds with 
theoretical evaluations, we start from the very beginning with a formalism which applies 
directly to this favorite particle candidate for cold dark matter. 
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The neutralino is defined as the hnear superposition 



(1) 



of lowest mass. Here 7, Z are the photino and zino states and H^, if^ ^i-re the higgsino fields, 
supersymmetric partners of the Higgs fields H^, The theoretical framework used in this 
paper for the supersymmetric model will be presented in Sect. 3. 

We write the differential event rate for elastic neutralino-nucleus scattering as 

= Nt^ dvf{v)v-—{v,ER). 2 

dEn Jv^^Er) dEn 



Nt is the number density of the detector nuclei, p-^jm-^ is the local (solar neighbour- 
hood) number density of neutralinos, j[v) is the velocity distribution of neutralinos in 
the Galaxy (assumed to be Maxwellian in the Galactic rest frame) evaluated in the 
Earth's rest frame, dajdE^ is the elastic differential neutralino-nucleus cross section, 
Eji = ml ^^v'^ {1 — cos 6*) / rriN is the recoil energy {6* is the scattering angle in the neutralino- 
nucleus center of mass frame), m^ed is the neutralino-nucleus reduced mass and mjv is the 
nuclear mass. Eq.(^ is written for a monoatomic material; its generalization to more com- 
plex materials is straightforward. 



A. Astrophysical and cosmological parameters 



Many astrophysical quantities, affected by large uncertainties, enter in the evaluation of 
the differential event rate. This is the case for the neutralino r.m.s. velocity and for its escape 
velocity, whose typical values are: f r.m.s. = 270 ± 24 km ■ s^^ [|r^, v^sc = 650 ± 200 km • s~^ 
pO| and for the velocity of the Sun around the galactic centre (f© = 232 ± 20 km • s~^ [@)- 



Also a large uncertainty concerns the value of the local dark matter density pi. A recent 
determination of p;, based on a flattened dark matter distribution and microlensing data, 
gives the range pi = O.Sllo^y GeV ■ cm~^ |21]. In particular the central value turns out 



to be significantly larger than the one previously determined, for instance, in Ref. ||22|| : 
Pi = 0.3 ± 0.1 GeV - cm~^. Furthermore, for any specific value for the local density of 
the total dark matter pi, one has to assign a value to the neutralino local density p^. To 
determine the value of p^ to be used in Eq.(||), we adopt the following rescaling recipe 
for each point of the parameter space, we take into account the relevant value of 
the cosmological neutralino relic density. When ^^^"^ larger than a minimal {Qh'^)min, 
compatible with observational data and with large-scale structure calculations, we simply 
put p^ = pi- When ^^^"^ turns out to be less than {Qh'^)fnin, and then the neutralino may 
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TABLE I. Values of the astrophysical and cosmological parameters relevant to direct detection 
rates. Vr.m.s. denotes the root mean square velocity of the neutralino Maxwellian velocity distribu- 
tion in the halo; Vesc is the neutralino escape velocity and V0 is the velocity of the Sun around the 
galactic centre; pioc denotes the local dark matter density and {Q,h?)min the minimal value of O/i^. 
The values of set I are the median values of the various parameters, the values of set II are the 
extreme values of the parameters which, within the physical ranges, provide the lowest estimates 
of the direct rates (once the supersymmetric parameters are fixed). 





Set I 


Set II 


Vr.m.s (km • S^^) 


270 


245 


yesc(km • s"^) 


650 


450 


^©(km-s-i) 


232 


212 


p;oc(GeV • cm-3) 


0.5 


0.2 


(rj/l )min 


0.03 


0.3 



only provide a fractional contribution Q^h"^ / {Qh'^)min = to Qh"^, we take = pi^. The 
value to be assigned to {Qh'^)min is somewhat arbitrary, in the range 0.03 < {Qh'^)min < 0.3. 

In Table I we give the two sets of astrophysical and cosmological parameters that we use 
in this letter. Set I corresponds essentially to the central values for the various parameters, 
whereas set II corresponds to those values of the parameters, which, within the relevant 
allowed ranges, provide the lowest detection rates (once the supersymmetric variables are 
fixed) . 



B. Cross-sections 



The differential cross-section -^^{v, Er) in Eq.(0) may conveniently be splitted into a 
coherent part and a spin-dependent one 



da I da \ I da \ 

(3) 

C \"'-^-K/ SD 



dEn \ dEn / „ \ dEpi 



When the WIMP is the neutralino, the coherent cross-section, 



..... ..^o. v^)c, originates from 

Higgs-exchange |^ and squark-exchange in the neutralino-nucleus scattering; the spin- 

•■dER' 



dependent cross-section {-^^)sd is due to Z-exchange and to squark-exchange The 
links among supersymmetric and nuclear degrees of freedom have quite different features in 
the coherent and in the spin-dependent cross-sections. Thus, let us consider the two cases 
separately. 
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1. Coherent x^nucleus cross-section 



In the case of the coherent cross-section, supersymmetric and nuclear degrees of freedom 

da ' 



decouple, thus {-j§-)c may be written as 



where is the total coherent cross- sect ion, conventionally integrated from zero up to E]^°'^ 
= 2ml^^v^ /rriN, q"^ = = 2m]yEji is the squared three-momentum transfer, and Fc{q) is 
a nuclear form factor defined as the Fourier transform of the nuclear matter distribution. 
Except for very light nuclei, which have a form factor Fc of exponential or gaussian type, a 
universal parametrization for Fc{q) is provided by 



gro 



where ji{qro) is the spherical Bessel function of index 1, s ~ 1 fm is the thickness parameter 
for the the nuclear surface, tq = (i?^ — Ss^)^/^ and R = 1.2 A^^^ fm (A is the atomic number). 



The cross-section o"^ may be written as 



a° = ^M|C^mLA^ (6) 

71 



where the quantity ( depends on the x ~ quark couplings mediated by Higgs particles and 
squarks . The full expression for ( together with the values used here for the relevant 



parameters are reported in Ref. ||2^. It is worth recalling that the simple dependence of 
(7(7 on A'^ is due to the equality (within a very good accuracy) of the ^-neutron and the 
X-proton couplings. 

From the factorization property inherent in the structure of Eq.(P) it follows that 
may be immediately converted into an equivalent x~ nucleon scalar cross section 

(nucieon) _ 1 + rn^/ mj^ 
'^'""'"'^ ~ 1 + ml/ml A^ 



where rup denotes the proton mass. Eqs.(H-0) enable a straightforward conversion of any 
upper limit on into an upper bound on cr^J^aiar''"^^ ■ This procedure allows a direct com- 
parison among the results of various experiments. It is obvious that this method applies 
also to any other WIMP which couples with equal strength to neutrons and to protons. 
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In Fig.l we present a summary of the current upper bounds, obtained by the method 
previously described. All experimental data have been treated with the same algorithm and 
with the same set of astrophysical parameters (set I above; {VLh?)min is not relevant here). 
We notice that the upper limit derived from the DAMA/Nal experiment provides the most 
stringent upper bound on the scalar neutralino-nucleon cross-section with an improvement 
of a factor 3-5 over the limits presently placed by Ge-detectors. The only exception is 
represented by a very narrow range in the neutralino mass around ^ 15 GeV. However, 
except for peculiar values of the other supersymmetric parameters, this range in is below 
the current lower bound on the neutralino mass [^,^ . 



2. Spin- dependent x^nucleus cross-section 



Let us turn now to the cross-section due to spin-dependent interactions. Using standard 

da 

dER 

da \ 16GI S{q) 



notations, we write the differential cross-section jp- for neutralino-nucleus scattering as 



SD 



where J is the nuclear spin and S{q) is given by: 

S{q) = alSooiq) + alSu{q) + aoaiSoi{q) (9) 



The functions Sij, = 0, 1) are rather complicated functions which depend on the nu- 
clear spin structure and on the neutralino composition; their definitions may be found in Ref. 



| 31| . The coefficients and ai are the isoscalar and isovector projections of the neutralino- 
neutron and neutralino-proton couplings dp, respectively, i.e. ao=ap+a„, ai=ap-a„. The 
full expressions for reported, for instance, in Refs. [p8| , p2| . In the case of the 

X-nucleon interaction mediated by Z-exchange one has 

ap,n= —{Ni - N!){ Y: T,,^<l)p,n (10) 

q=u,d,s 



where T^q denotes the third component of the quark weak isospin, Ag is the fractional spin 
carried by the quark q in the appropriate nucleon and iVs, N/^ are the Higgsino coefficients 
in the expression for x (see Eq.(P). 



The cross-section of Eq. (§) may be conveniently rewritten as |^ 
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da \ 16G% .9^/^ ^S(q) _ . 

> - ^-mjvAV(J + l)-)^, (11) 



where the static value of S{0) has been expressed as 

2J+1 .2 



^(0) = — AV(J + 1) (12) 
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with 



^ ^ ttp < Sp > +a„ < Sn > (13) 



< Sp > and < Sn > are the contributions of the total proton and neutron spins in the 
nucleus. Another commonly used notation is Ap(„) = A/ap(„). 

It is important to recall that both A and S{q)/S{0) depend on nuclear properties as well 
as on the neutralino composition. Thus, in spin-dependent interactions nuclear degrees of 
freedom and supersymmetric degrees of freedom do not decouple in the most general case. 
This feature does not allow an extraction of a universal {i.e. independent of the neutralino 
composition) spin-dependent nucleon cross-section in the same way as was done in the 
coherent case. 

An exception is provided by spin-dependent interactions mediated by Z-exchange only. 
This property is explicitely exhibited by the expressions of Eq.(p!OD for ap, a^, which imply 

K/ap]z=(Y.T3,Aq] /(j:T,,Aq] (14) 



In order to compare experimental data due to various detectors in a way similar to the 
one previously employed for the coherent interactions, in the rest of this section we expand 
the formalism in the case of a pure Z-exchange contribution to the spin-dependent scattering. 
Under this assumption, the spin-dependent cross-section may be cast into the form 



where Xp^z is defined as 
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A 



J 



<Sp> +— <Sn> 



(16) 



an/ap = [a„/ap]z 



and 



SD 



' S{q) 
S{0) 



(17) 



a.n/ap=la„/ap]z 



Thus we may further define a total point-hke spin- dependent cross-section (convention- 
ally integrated from zero up to (Ejijmax) as 

r TP \ ' ' ■'- 

^SD ~ [ERjmax 



^^'-mUN! - Nlf (E^3,Ag) Aj^J(J + 1), (18) 



TT 



in terms of which we may obtain an equivalent proton spin-dependent cross-section 

[proton) _ 1 + m'^/mlf 3 CTgji , . 

""^^ " l + ml/mlAXl^JiJ+l) ^ ^ 



Therefore, under the assumption of a pure (or dominant) Z-exchange, Eqs.( P^P^ allow 
the conversion of any experimental upper limit on into an upper bound on ag^^"^\ in 
much the same way as for the coherent case. However, still an important difference between 
the two cases concerns the features of the nuclear properties. Whereas in the coherent case 
these are simply taken into account by the form factor Fc of the universal analytic form 
of Eq.(H), in the spin-dependent one the nuclear properties have to be calculated for each 
nucleus with an appropriate nuclear model. In our analysis of the spin-dependent cross- 
sections we have employed the most recent nuclear-physics calculations for the structure 
functions Sij{q) for the nuclei of interest. From the previous formulae it is clear that from 
these Sij's both the static value A and the form factor Fsd are readily derived. 

The nuclei with non-vanishing spin, relevant for our analysis, are ^^Na, ^^Ge, ^^''I, ^^^Xe. 
For these nuclei we give in Table II the values for J{J + 1) used in our analysis P,|3|,|36 



compared with some of the previous determinations pi|P5[ . We have set a„/ap = [an/ap]z, 
and we have assigned to this ratio the value [an/ap]z = —0.85, which is obtained from 
Eg . (p!0D using for the Ag's the estimates of Ref. . 
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TABLE II. Values of the nuclear spin-dependent static properties. In the first column are listed 
the nuclei with non- vanishing spin considered in the present paper. Their nuclear spins are given in 
the second column. In the third, fourth and fifth columns are listed the values of < Sp > , < Sn > 
and Xp^J{J + 1) derived from the following references: Ref. Q for ^^Na and ^^^Xe, Ref. |36] for 
''^Ge, and Ref. ||2[^ for ^^"^I. In the sixth column is given the ratio {\{n),z/ ■^og)'^- the numerator 
stands for the A factor for proton or for neutron, according to whether protons or neutrons are the 
odd species in the relevant nucleus {Xp^z for ^^Na and ^^''l, and An,z for ^^Ge and ^^^Xe), and its 
value is obtained from the Refs. P,P,|36[ ; the denominator stands for previous evaluations with the 
odd-group approximation for ^^Na, ^^Ge and ^-^^Xe [34| and with the interacting boson fermion 
model for ^^''l (B) and (N) denote that the nuclear-physics calculations of Refs. |2|,^ were 

performed with the Bonn A or with the Nijmegen II nucleon-nucleon potential, respectively. 



nucleus 


J 


<Sp> 


< Sn > 


K,zJiJ + '^) 


'^l{n),z/'^OG 


23Na 


3/2 


0.248 


0.020 


0.089 


2.2 


73Ge 


9/2 


0.030 


0.378 


0.103 


2.2 


1271 (B) 


5/2 


0.309 


0.075 


0.084 


3.6 


1271 (N) 


5/2 


0.354 


0.064 


0.126 


5.5 


i29Xe (B) 


1/2 


0.028 


0.358 


0.229 


2.5 


i29Xe (N) 


1/2 


0.013 


0.300 


0.177 


2.0 



It is noticeable that the new evaluations of the static properties provide larger values for 
A as compared to the previous ones. We remark that these quantities are somewhat sensitive 
to the type of the nucleon-nucleon interactions used in the calculations. 

The form factor Fsd for the various nuclei has been calculated using the nuclear-physics 
results of Refs. 0J^,^. The relevance of these finite-size effects depend on the range of 
the variable qR which is employed for the extraction of the upper bounds on the differ- 
ential counting rates. This, in turn, depends (apart from the nuclear size) on the other 
experimental features: energy threshold and level of the counting rates. For "^^Na. and 
^^Ge finite size effects are negligeable, since the relevant qR values are 0.27 < qR < 0.55 
for ^^Na and 0.75 < qR < 1.4 for '^^Ge. On the contrary, these effects have to be 
properly taken into account for ^"^"^l and ^^SXe, whose qR ranges are: 2.2 < qR < 4.4 
for i27j and 1.6 < qR < 3.5 for ^^SXe. In Fig. 2 are shown the plots of the function 
[S{q)/ap]z = [S{q)/ap]a„/ap=ian/ap]z for these two nuclei, calculated from the structure func- 
tions Sij's of Refs. 0,^ for the two different nucleon-nucleon interactions previously men- 
tioned. Again we have set an/ dp = [an/cip\z = —0.85. It is worth noticing two main features: 
i) the shape of S{q) is not very sensitive to the form of the nucleon-nucleon interaction, ex- 
cept for i29Xe for qR > 4; a region which is, however, beyond the qR-range employed in our 
extraction of the upper limits, ii) the function S{q) for ^^^j shows a plateau for qR > 3, then 
in a region which is important in our analysis. 



When we apply the procedure previously outlined {i.e. Eqs.(|T5|-|T9|)) to the counting 
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rates of the experiments of Refs. |lT|,[5|j9|- [TT| , p!3| we finally obtain the results shown in Fig.3. 
For the quenching factor Q the following values have been used: Q = 0.30 for ^^Na [[li|7|,pj]| , 
0.25 for 73Ge [|, Q = 0.09 for f|Il,|7|JTl and Q = 0.65 for ^^^Xe 



12,14 . From this 



Q = 0.25 for ^^Ge [|, Q = 0.09 for ^"^n [|l],0JTl| and Q 
figure it turns out that the data of the DAMA/Nal provide the most stringent upper bound 
also for the spin-dependent cross-section a^^"*""''. It is noticeable that, because of the new 
higher value for A and of the property ii) discussed above, the constraint from the I-nucleus 
turns out to be significantly stronger as compared to results of previous analyses. 



III. EVALUATION OF CROSS-SECTIONS IN THE MSSM 



The supersymmetric model employed in this paper is the so-called Minimal Supersym- 
metric extension of the Standard Model (MSSM) |3^, which provides the most convenient 
phenomenological framework at the electroweak scale (Mz) and does not assume too strong, 
arbitrary theoretical hypotheses (for instance, universality conditions for the soft scalar 
masses at the Grand Unified scale Mqut)- Various properties (relic abundances and detec- 
tion rates) of relic neutralinos have been analyzed in the MSSM by a number of authors. 
Some of the most recent ones are given in Refs. |50J5^H0-H2 . 



The MSSM is based on the same gauge group as the Standard Model and contains the 
supersymmetric extension of its particle content. The Higgs sector is modified as compared 
to that of the Standard Model, because it requires two Higgs doublets Hi and H2 in order 
to give mass both to down- and up-type quarks and to cancel anomalies. Thus the MSSM 
contains three neutral Higgs fields: a pseudoscalar and two scalar fields. The Higgs sector is 
specified at the tree level by two independent parameters: the mass of a physical Higgs field 
(we use the mass of the neutral pseudoscalar boson) and the ratio of the two vacuum 
expectation values, usually defined as tan/5 = V2/V1 =< H2 > / < Hi >. Other important 
ingredients of the model are, apart from the Yang-Mills Lagrangian, the superpotential, 
which contains all the Yukawa interactions and the Higgs-mixing term fiHiH2, and the 
soft-breaking Lagrangian, containing the trilinear and bilinear breaking parameters and the 
soft gaugino and scalar masses. 

The theoretical model contains an exceedingly large number of parameters, unless some 
assumptions are introduced. It is customary to introduce the following conditions (at the 
Mz scale) : i) all trilinear parameters are zero except those of the third family, which however 
have a common value A, ii) all squarks are taken as degenerate: rrig. = tuq, except for m^, 
iii) all sleptons are taken as degenerate with m^- = mo, iv) the gaugino masses are assumed 
to unify at Mqut] thus in particular Mi = (5/3) tan^6'i^M2. 

After these conditions are applied, the number of independent parameters reduces to 
seven. We choose them to be: M2, /i, tan/3, m^, mo, A. In terms of these parameters, 
the supersymmetric space is further constrainted by the following requirements: 1) all ex- 
perimental bounds on Higgs, neutralino, chargino and sfermion masses are satisfied (taking 
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into account also the new data from LEP2 p9| , ^ , ^ ) , 2) the neutrahno is the Lightest Su- 



persymmetric Particle (LSP) (i.e., regions where gluino or squarks or sleptons are LSP are 
excluded), 3) constraints due to the 6 — > s + 7 process are satisfied, 4) the neutralino 
relic abundance does not exceed the cosmological bound, i.e. fl^h'^ < 1. 

The results presented in this paper for the evaluations of neutralino-nucleon cross sections 
within MSSM have been obtained by varying the independent parameters of the model in 
the following ranges: 10 GeV < M2 < 500 GeV, 10 GeV < < 500 GeV, 60 GeV < 
tua < 500 GeV, 100 GeV < mo < 500 GeV, 100 GeV < m^- < 500 GeV, -3 < A < 
+3, 1.1 < tan/5 < 50, and taking into account the experimental constraints previously 
mentioned. The ensuing explored range for is from its lower bound up to 250 GeV. 

Within the present theoretical framework we have evaluated the cross-sections o'l^atar"^'^ 
and 0"^^°*°"^ The results for <J^^aial°^\ multiplied by the scaling factor are shown in 
Fig. 4. Allowed neutralino configurations stay inside the outer contour line. To provide more 
information, also the locations of configurations with the representative values: 60 GeV < 
< 70 GeV and tan /5 = 10, 50 are displayed. The thick curve represents the experimental 
upper bound to cy'^^aiil°^\ deduced by using the set I of Table I; thus this curve is just a part 
of the one in Fig.l and is due to the new data from the DAMA/Nal experiment. This figure 
shows that a large number of neutralino configurations are above the current experimental 
bound. It follows that the present experimental sensitivities have reached a level of great 
significance for the possibility of investigating realistic supersymmetric models. However, 
this result cannot be used to disallow those configurations which provide a ■ cricatar"^^ above 
the experimental limit, because of the uncertainties affecting the astrophysical parameters. 
If we assign to these parameters their most conservative values, i.e. the values of set II of 
Table I, we obtain the results displayed in Fig. 5. This Figure shows that, even in the very 
conservative scenario represented by set II, for the first time the boundary of the physical 
region for neutralinos is reached by direct detection measurements. From our results it 
turns out that the neutralino configurations currently accessible by direct detection are 
those pertaining to light Higgs masses: rriA — 60 GeV and large tan/3: tan/3 ~ 40 — 50. 
These same configurations are at present under investigation at LEP2 (at y/s = 172 GeV). 

A similar analysis, performed for the spin-dependent case, shows that the theoretical 
predictions within MSSM for a^^"*""^ are still much below the current upper bound of Fig. 3 
(by more than two orders of magnitude). This is expected, in view of the much reduced 
strength of spin-dependent effects versus the coherent ones. 

A few comments are in order here. Investigations of fundamental particle-physics prop- 
erties by non-accelerator experiments usually suffer from a number of serious uncertainties. 
This is the case for the problem under discussion, as we have already recalled above, in 
connection with the astrophysical and cosmological parameters. For this reason, in our 
conclusions we have also conservatively taken into account the set II for these parameters. 
However, it is remarkable that, even in this unfavorable scenario, the present sensitivity of 
experiments for direct WIMP search turns out to be at a level of great physical interest for 
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supersymmetric relic particles. Significant improvements in the determination of fundamen- 
tal cosmological parameters (such as the average matter density and the Hubble constant) 
are expected in the near future; these will allow to draw stronger conclusions from investi- 
gations similar to the one presented here. It is also to be remarked that further substantial 
improvements in the sensitivities of experiments searching for WIMPs are expected, thus 
the perspectives for this kind of experimental investigation are very encouraging. 



IV. CONCLUSIONS 



Using the results of an experiment for direct detection of Weakly Interacting Massive 
Particle (WIMP), which employs a large- mass, low-background Nal(Tl) detector at the Gran 
Sasso Laboratory, DAMA/Nal and new nuclear-physics evaluations of spin-dependent 
structure functions we have shown that these new experimental data provide the most 
stringent upper bounds both to scalar and to spin-dependent WIMP-nucleon cross-sections 
(except for a very narrow range for WIMP mass around 15 GeV). Applying these results to 
the most favorite candidate for cold dark matter, the neutralino, we have proved that, for 
the first time, data from an experiment of direct search for WIMPs allow the exploration of 
a region of the supersymmetric parameter space not yet fully investigated at accelerators. 
This result has been obtained by using for the astrophysical and cosmological parameters a 
very conservative set of values. We may also conclude that the experiment of Ref. ^ can 
allow an investigation of modulation effects with significant sensitivity. 
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Figure Captions 



Figure 1 - Upper bounds on o'l'^aiar"'^^ ■ The dot-dashed hne denotes the hmit obtained 
by combining data of Ge-detectors; the branch for < 10 GeV is from results of Ref. |]10 



the one for 10 GeV< < 20 GeV is from results of Ref. P and the one for > 20 GeV 
is from data of Ref. |Q. The solid and the long-dashed lines denote the bounds from the 
results of Nal-detectors: the solid line is from data of Ref. |]l|, the long-dashed line is from 
data of Ref. |TT|. Dotted and short-dashed lines are from data of a Xe-detector []T3[ and of 



a Te02-detector [|15[, respectively. 

Figure 2 - Structure function [S{q)/ap]z, obtained from the results of Refs. versus 
qR. The two lower curves refer to ^^^Xe with the Bonn A (long-dashed line) and with the 
Nijmegen II (dot-dashed line) nucleon-nucleon potentials. The two upper curves refer to ^^''I 
with the Bonn A (short-dashed line) and with the Nijmegen II (solid line) nucleon-nucleon 
potentials. 

Figure 3 - Upper bounds on The solid and the dashed lines denote the bounds 

from the results of Nal-detectors: the solid line is from data of Ref. [|l|, the dashed line is 
from data of Ref. ]ll|]. The dotted line is from data of a Xe-detector |jl3| and the dot-dashed 



line is from Ge-detectors P,PJTO 

Figure 4 0"^"^;^'^°"^ multiplied by the scaling factor ^, versus the neutralino mass. Al- 
lowed neutralino configurations stay inside the outer contour line. Configurations with 60 
GeV < tua < 70 GeV and tan (3 = 50 (tan P = 10) stay in the up-left to bottom-right 
hatched region (in the doubly hatched region). The thick line denotes the upper bound 
obtained from data of Ref. |I| . The set I of Table I for the astrophysical and cosmological 
parameters has been used. 

Figure 5 The same as in Fig.4 except that set II of Table I for the astrophysical and 
cosmological parameters has been used. 
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